Analytical formula connecting entangled state and the closest disentangled state 
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The separable state closest to a given entangled state in the relative entropy measure is called 
the closest disentangled state. We provide an analytical formula connecting the entangled state and 
the closest disentangled state in two qubits. Using this formula, when any disentangled state (a) 
located at the entangle-disentangle boundary is given, entangled states to which a is closest can be 
obtained analytically. Further, this formula naturally defines the direction normal to the boundary 
surface. The direction is uniquely determined by a in almost all cases. 

PACS numbers: 03.67.-a, 03.65.Ud 



Quantum entanglement is the most striking feature of 
quantum mechanics. Several measures have been pro- 
posed to quantify the amount of entanglement . One 
of the important measures is relative entropy of entan- 
glement 2, 3], which is defined as 

Er(q) = min S(g\\o~) = min [Trying — Trainer] , (1) 

where the minimization is performed over all density ma- 
trices in the set of disentangled states (In denotes the 
natural logarithm throughout this paper). The states 
achieving the minimum are called the closest disentan- 
gled states. The efficient algorithm for the numerical 
minimization has been proposed in [4j . However, the an- 
alytical minimization is difficult in general, and the an- 
alytical solutions only for several classes of states have 
been shown so far. 

In this paper, we provide an analytical formula con- 
necting entangled states and the closest disentangled 
state in two qubits. Using this formula, when any dis- 
entangled state (cr) located at the entangle-disentangle 
boundary is given, entangled states g to which a is clos- 
est can be obtained analytically. Further, this formula 
naturally defines the direction normal to the boundary 
surface when relative entropy is regarded as the distant 
measure. The direction is uniquely determined by a in 
almost all cases. The discussions in this paper are limited 
to the case of two qubits. 

The strategy of the minimization is based on the fol- 
lowing fact: Any full rank density matrix can be trans- 
formed into a unique Bell diagonal state by applying a 
suitable local filtering operation 0, |(| . The explicit for- 
mulation using the Lorentz matrices has been shown in 
, and that using the Wootters basis has been shown 
in 0. The entanglement manipulation by local filter- 
ing on a single copy was originally considered in 0, . 
Further, the connection between the Wootters basis and 
Lewenstein-Sanpera decomposition |lCj has been shown 
in 0| • According to these results, any full rank density 
matrix can be written as q—j^(Fa ® Fb)qbd(f\ ® -^s)) 
where N — Ti[(Fa ® Fb)qbt> (F\ <g) F B )} and QBD = 

~Y^i=oPi\ e i)( e i\ i s the Bell diagonal state. £>bd c & n be 
always chosen in a canonical form, in which pg is max- 
imum (po > pi,P2,Pa) and \e%) is a fixed set of Bell 
basis, since suitable local unitary operations transform 



any Bell diagonal state into the canonical form 0, ^| . 
On the other hand, g can be always expressed using 
the Wootters basis as g = ^2 i=0 W<j>i){4>i\ 0, where 
(4>i\4>j) — ((f>i\(f)j) = 5ij defining the tilde operation as 

= ((72 ®o-2)\ip*)- The Bell basis is one of the Wootters 
basis by itself (adding a suitable global phase if neces- 
sary), and hence |ej) = |ej). The concurrence of g is 
C = A — Ai — A 2 — A 3 . Although the actual definition of 
the concurrence is max(Ao — Ai — A2 — A3, 0) , we define 
it here such that C can take negative when g is disentan- 
gled. Further, we can put det Fa = det Fb = 1 without 
loss of generality. In this choice, \<j>i) = {Fa ® Fb)^) 
(using A T a2A= (det A)o-2, it can be checked that \4>i) = 
{F A ®F B )\ei) satisfy (<tn\fa) = 5 ij ) and A i =p i /7V [Ell 

Considering these facts carefully, the following struc- 
ture of the Hilbert space is revealed: The total 15- 
dimensional space is divided into subspaces in such a 
way that only the states in the same subspace can be 
transformed to each other by local filtering (we only 
consider the full rank density matrices for the sake of 
simplicity [EH)- Each subspace contains the canon- 
ical Bell diagonal states having the same pi only (we 
do not care about the exchange among p\, P2, and ps). 
Therefore, the subspaces can be specified by three pa- 
rameters (pi with ^2 Pi — 1), and each subspace consti- 
tutes a 12-dimensional manifold. Further, since C > 
for entangled states and C < for disentangled states, 
and the local filtering does not change the sign of C (see 
Theorem 1 in 0), we can define boundary subspaces, in 
which all the states have just zero of C. The bound- 
ary subspaces are specified by two parameters, and the 
assemble of the boundary subspaces constitutes the 14- 
dimensional boundary surface separating the region of 
entangled states and disentangled states. 

Then, we perform the minimization of S , (f||cr) sepa- 
rately: First we obtain the extremal conditions by mini- 
mizing within a fixed subspace, and second we obtain the 
conditions to determine the boundary subspace contain- 
ing the closest disentangled state. Finally we solve the 
equations obtained. Hereafter, the entangled state and 
the corresponding closest disentangled state is denoted 
by g and a, respectively. Further, we first assume that 
a is full rank. The case of a with lower rank will be 
discussed in the last part of this paper. 
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We have already performed the first step and obtained 
the set of equations a must satisfy Q. We briefly 
repeat the results here. The states obtained from a 
by local filtering belong to the same subspace as a by 
its definition. Let us consider Bob's local filtering of 
a' = ® e m ' 5 l 2 )a(I ® e m ' 5 1 2 ) , where t is a real param- 
eter and a ~ (a%, 02, 03) is the vector of Pauli matrices. 
In order that 5 l (g||(T) is minimum, the linear coefficient 
of t in the expansion of Trplncr' must be zero. As the 

result, using ln(X + tY) = \nX + 1 f^x^Y^dz H 

we obtained 

Tr(7 ® a t )Za + TrZ(I ® <t 1 )<t = fori = 1,2, 3 (2) 
where 

r°° 1 1 

Z= / q dz-I (3) 

J a + z a + z 

(Eq. (J2J) holds for i = also, but we omitted it since 
TrZa = is obvious from Eq. 10 )• In the same manner, 
considering Bob's local unitary transformation of a' = 
(I <g> e itn -^ 2 )a(I ® e - it ™- <? / 2 ), we obtained 

Tr(J ® <ti)Zo - TrZ(I ® a^a = fori = 1,2, 3. (4) 

From Eqs. (0) and 10} , we have 

Tr(J 0^)^ = 0, (5) 

which is sufficient for TiZ(I ® o'i) (T = since Z is Her- 
mitian. Therefore, a must satisfy Eqs. © and Alice's 
counterparts 

Tr(cri ® 7)Zcr = 0. (6) 

These are fundamental equations for the closest disen- 
tangled state (and even for the closest positive partial 
transposed states), because these equations must hold 
in any multi-party systems and any higher dimensional 
systems [l4j. 

The second step is to obtain the conditions to deter- 
mine the boundary subspace containing a. The canonical 
Bell diagonal state in the boundary subspace must have 
zero of C = po — Pi — P2 — P3, but ^2 Pi = l due to the 
normalization. As the result, po must be 1/2 (this is also 
obvious from the fact that the concurrence of the Bell 
diagonal state is 2pg -10 El)- Therefore, cr, which 
belongs to the boundary subspace, must be obtained by 

local filtering from <t B d = ||eo)(eo| +2i=iP»l e i)( e »l> and 
hence a must be written as 

a = -^{F A ®F B )a B D{Fi®Fi) 

= ^I0o)<0oi+E|i^)^i' ( ? ) 

i=l 

where N is again normalization, and the set of |</>i)'s is 
the Wootters basis of a. 

Adopting this parameterization, to determine the sub- 
space containing a is equivalent to determine p\, P2, 
and j»3 under the constraint of p\ + P2 + P3 — 1/2 and 



< Pi < 1/2. However, we assumed that a is full rank, 
which ensures that the minimum of i5(g||o") is achieved 
at 0<Pi < 1/2. Therefore, we can simply minimize 

f = -Trp\na + l( Pl +p 2 +p 3 -l/2), (8) 

with I being a Lagrange multiplier. Considering the small 
change of pi (note that N also contains pi), we obtain 

^TiZ\^){4n\=l for i = 1,2, 3. (9) 

Further, multiplying pi on both side of the above equa- 
tions and adding them, we find 

£^TrZ|&)<&| =Tr^-^Tr^ )^ | = L, (10) 

i=l 

and we finally obtain the extremal conditions as 

(4> i \Z\<l> i ) + ((t> \Z\<t> )=0 fori = 1,2, 3. (11) 

There are three equations in the above, in spite that the 
boundary subspaces are specified by two parameters as 
mentioned before. This is because Eqs. Ijllfl includes the 
trace condition of Trp = Trcx. 

Then, all equations of the extremal conditions are Eqs. 
©, ®, and |TT} with pi + p 2 + Pi = 1/2 (and thus 
Ao — Ai — A2 — A3 = 0). The total number of these equa- 
tions is 16 including the normalization condition. These 
must be solved for any given g in order to obtain a true 
closed formula for the closest disentangled state (and for 
Fr{q))- However, since it is a hard problem unfortu- 
nately, we follow the strategy in [l^. Namely, for any 
given cr, we obtain the set of g for which 5(gi||cr) is min- 
imum. We can obtain it in an analytical form as shown 
below. 

Let us assume that a with C — is given. Here, we still 
assume that a is full rank. It is convenient to express a 
in two representations: 

a = ^A ii |^)(^|=2r ii |i)01, (12) 

ij ij 

where |<fo)'s and |i)'s are the Wootters basis and eigen- 
states of cr, respectively, and hence Ay = XiSij and 
Tij = li&ij with 7j's being eigenvalues of a. These 
two basis sets are connected through a unitary matrix 
as ^/Xi\4>i) = Yjj V i*j\/lj\j) 0- Since l<^)' s are non- 
orthogonal in general, it is also convenient to introduce 
the matrices describing the non-orthogonality as follows: 

Qy = <&|&>, n« = (13) 

Using the completeness of I = J^i \<t>i){4>i\ = J2i \<t>i){4>i\ 
III and I = J2ij Qij\4>i) {h I = Ey % \<Pi) (&|, it is 
easy to check that Hermitian Q and n satisfy QQ T = I, 
n = Q*, and QII = T1Q = L Further, 

r = v^VaqVav = Uy/Qky/Qu\ (w) 
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where we introduced the unitary U by the singular value 
decomposition of U\[Q yAV — y/T. Then, we shall obtain 
g in the same representation as a. Namely, using \<f>i) and 
\i) of a, 



(15) 



Hereafter, the coefficient matrix in the Wootters ba- 
sis and eigenstates is denoted by suffix A and E, re- 
spectively. Two representations are connected through 

Viewing the form of Eqs. (JSJl, ©, and (|ll|l . it is found 
that, when g is a solution, g' = xg + (1 — x)a is also 
a solution, since Z — ► xZ for q — * xg + (1 — x)a. This 
implies that the normalization of Z can be taken as a 
free parameter, and we express Z as 



Zj) = / — 

Jo it ^ 



r>E 

+ z ij 



1 



Then, we have xW^ : 



7j 



(16) 



G] — /, where AoB is the 



Hadamard product defined as [A o Z?]y = AijBij, Gij — 
\/Gij, and 



Gy = 



for 7i=7j 
in - .i,,-, for 7 ^ 7j . 



7» 

7i~7j 



(17) 



By inverting the above (we can safely invert since Gy > 
by 7i > 0), we have 



i? B = r + i[w B oG], 



(18) 



and therefore the problem is reduced to find W E satisfy- 
ing the set of all extremal conditions. 

For this purpose, it is convenient to express Z in the 
\cf)i) basis as Z — ^2xWij\^>i)(<j>j\. Then, the extremal 
conditions of Eqs. Ijllfl become quite simple: 

Wii=-W 00 fori =1,2,3. (19) 

Since the normalization of W was absorbed in x, we 
adopt Woo = 1 and the diagonal elements of W are de- 
termined to be {1,-1,-1,-1}. The remaining extremal 
conditions are Eqs. 10) and (J5J), those are 



lm 

E< 

l in 



)W ml X l = 
)W ml X l = 0. 



(20) 



On the other hand, by virtue of F\Pa = F\oiF* A oi = / 
n [l3| (we put det F A = det F B = 1) , we find 



(<h\I ® tn\<h) = {ei\F\F A ®Fla i F B \ei) 
= Ti(I ® FlaiF B )\ei){ei\ 

= ^ra i <T 2 F^a 2 F t B = 2 Trcr » = °i 

and therefore, Eqs. I|20|l are the set of linear equations 
only for the off-diagonal elements of W. As the result, 
W m ^i=0 obviously satisfy Eqs. (|2Tj|l . After all, 



W = diag{l, -1,-1,-1} 



(21) 



is a solution of all extremal conditions. Since W E = 
U \/TlWVTlU' < , we finally obtain a relation between g and 
a in an analytical form as 



R E 
A E 



r + xA E _ 
[(uVnwVmi^) o g] 



(22) 



with G by Eq. (JT^, W by Eq. (EH, and n by Eq. fT3). 

It should be noted that the phases in U, as well as V, 
should be chosen so that {<f>i\<f)j) =Sij. In order to avoid 
the complication, it is convenient to express it in the 
Wootters basis, that is 



R = A 



a a = Vuu ii [(uVTiwVuu'') o g]uVu. (23) 

In this expression, any U satisfying U y/QA^/QU^ = T 
can be used safely. From the relation, it can be seen 
that the non-orthogonality of the Wootters basis plays 
an important role. 

The trace of the matrix 6 = X^jj A^|z)(j| is confirmed 
to be zero as 

TtA e = ^[uVnwVTiu%a i = TiuVTiwVTiu^T 

= TrVFA = A - Ax - A 2 - A 3 = 0, (24) 

and S obviously plays the role of the direction normal to 
the boundary surface at the position of a, when S(g\ \a) is 
regarded as the distant measure between states. Which 
side does the entangled region spread? For an infinitesi- 
mally small by treating A A as the perturbation, the 
diagonal elements of A A contribute to the concurrence of 
g in the lowest order of and then 



C(g) ~ 
AG = 



xAG 



Tr A A W = V | [UVUWVUU^] i . 1 2 G tJ > 0, 



since Gij > and U VTlWVTiU^ is not zero operator 
(it is invertible). Therefore, the entangled region always 
spreads in x>0. When a is full rank, we can always find 
the positive g for some infinitesimally small \x\. This 
implies that there are both entangled states and disen- 
tangled states in the vicinity of a. Therefore, the full 
rank a with G = is the true boundary state separating 
the entangled and disentangled region. Further, Er(q) 
for |x|<Cl is also calculated as 

S{g\\a) - xTy(Z + 1)5 = x 2 TtWA a = x 2 AC. (25) 

It is important to note that there are entangled states 
satisfying Eqs. iJSJ), (JBJ), and (|llfl for which a is not opti- 
mal. However, the solution we obtained are not the case: 
cr is indeed optimal for g = a + xS with any x > 0. This 
can be proven as follows. The equations we solved take 
into account any small deviation from a along the bound- 
ary surface, and the solutions are ensured to be extremal 
along any path on the boundary. This extremeness does 
not immediately imply the global minimum of S^pHc) 
among the all disentangled states by itself. However, for 
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our solution, S(g\\a + xS) > S(g\\a) holds for x < (in 
the vicinity of er), and hence the relative entropy inside 
the boundary is larger than that on the boundary. In 
this situation, considering the semi-continuity and con- 
vexity of the relative entropy |18|| , it can be shown that 
the extremeness along the boundary must be local min- 
imum among the disentangled states around er, and it 
must be global minimum |3(. Further, the assumption 
we first made that g has er of full rank is also justified in 
the same reason. 

It should be noted further about the uniqueness of W . 
(if it is not unique, there might exist entangled states 
other than those we obtained, for which er is optimal). 
Since W is Hcrmitian, the number of real parameters of 
the off-diagonal elements is 12, and the number of lin- 
ear Eqs. I|20|) is also 12. Therefore, Eq. I|21|l is a unique 
solution unless some of Eqs. lf2*U|l are dependent to each 
other. We have numerically generated random cr's and 
confirmed that Eqs. I|2(J|) are always independent. How- 
ever, such Monte Carlo type calculations cannot pick up 
the states of measure zero, such as those having some 
special symmetry, and the other solutions g have not 
been completely excluded. Further, it is an open ques- 
tion whether the solution we obtained is enough to ex- 
press whole entangled states (does a + \x\S span whole 
entangled states, when a moves on the boundary?) 

In so far, we solely discussed the case that er is full rank. 
In the case of er with a lower rank (denoted by <tl), let us 
consider a sequence constituted by all full rank er n 's con- 
verging as <r n — ><7£. We can always find such sequences 
because there always exist full rank boundary states in 
the vicinity of trz, (if of er^ is zero, adding infinitesi- 
mally small portion of \4>k){4>k\ to ax). In this sequence, 
g n obtained from Eq. (|22|l also converge as g n — > goc 
(since Gy is non-singular). Then, the continuity of the 
relative entropy of entanglement shown in |19| ensures 
that E^goa) = S{g ao \\uL), and ox is one of the closest 
disentangled states of g^. After all, Eq. (|2*2*)l must hold 



under such limiting sequences even for low rank ex. 

Finally, we show some simple examples of Eq. <|23ll . 
The Wootters basis can be easily constructed by \<j>i) = 
(Fa^Fe)^). For Fa = Fb=I, a becomes a Bell diagonal 
state, and [</>,) = |e^) constitute the orthogonal set. Then, 
Q = II = U = 1 and T = A, and we obtain 

A A = diag{l/2,-A 1 ,-A 2 ,-A 3 } 

R A = diagji^a-^Ax, (1-2^2,(1-^3}. 

For the choice of Fa = F b = diag{£, l/t}, |e ) —i\4> + ), and 
\ei) = (I®ai)\<f) + ), where |e6+) = (00) + |ll))/\/2, we have 



Q = n* = - 



/ t 4 + i/t 4 



\tt 4 -i/t 4 



-%t 4 + i/t 4 



> (26) 



t 4 + l/t 4 



and the limiting sequence of Ai = A2 
gives T = A diag{t 4 ,0,0, l/t 4 }, and 

A A cx diagjl, 0,0,-1} 



and A 3 



R A = diagjl+a/,0,0, l-x'}. 



Ao 



(27) 



These examples reproduce the correct relation between 
Bell diagonal g and er , and between pure state and its 
closest disentangled state of rank 2 Q , respectively. 

In summary, we provided an analytical formula con- 
necting er and g in two qubits [Eq. I|22|l ]. Using this 
formula, when a with C — is given, g for which Seller) 
is minimum can be obtained analytically. Further, this 
formula naturally defines the vector (5) normal to the 
cntangle-disentangle boundary surface. The normal vec- 
tor is uniquely determined by a in almost all cases. 
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